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ABSTRACT 

The exact nature of r\ Carinae is still an open issue. Strict periodicity in the light 
curves at several wavelengths seem to point out to a binary system, but the observed 
radial velocities, measured from space with high spatial resolution are in conflict with 
the ground based observations used to calculate the binary orbit. Also, the observed 
2-10 keV X-ray flux is much larger that what is expected from a single star, and 
favors the wind-wind collision hypothesis, characteristic of high mass binary systems. 
However, to explain the duration of the dip in the light curve by wind collisions, it is 
necessary to postulate a very large increase in the r\ Carinae mass loss rate. Finally, the 
optical and UV light curves are better explained by periodic shell-ejection events. In 
this paper we conciliate the two hypothesis. We still assume a binary system to explain 
the strong X-ray emission, but we also take into account that, near periastron and 
because of the highly eccentric orbit, the wind emerging from r\ Carinae accumulates 
behind the shock and can mimic a shell-like ejection event. For this process to be 
effective, at periastron the secondary star should be located between rj Carinae and 
the observer, solving also the discrepancy between the orbital parameters derived from 
ground and space based observations. We show that, as the secondary moves in its 
orbit, the shell cools down and the number of available stellar ionizing photons is not 
enough to maintain the shell temperature at its equilibrium value of about 7500 K. 
The central part of the shell remains cold and under these conditions grain formation 
and growth can take place in timescales of hours. We also calculated the neutral gas 
column density intercepting the line of sight at each point of the orbit near periastron, 
and were able to reproduce the form and duration of the X-ray light curve without 
any change in the r\ Carinae mass loss rate. This same column density can explain the 
observed Ha light curve observed during the 2003 event. 
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1 INTRODUCTION 

r\ Carinae is a peculiar star surrounded by a dense cloud of 
gas and dust, formed in several episodes of mass ejection, 
typical of LBV stars. Its mass, calculated from the bolo- 
metric luminosity, corresponds to approximately 120 Mq, 
making r\ Carinae, if single, the most massive star in the 
Galaxy. During the last years strong evidences of the system 
binary nature were reported: a 5.52 year periodicity in th e 
light curve at infrared wavelengths i Whitelo cj^^ilJ| W9^j . 
in the Hel A10830 A high excitation line \D aminelil 1 1 9961) 
and in the radial velocity of the Pa7 line, the latter pointing 
out to a highly eccen tric orbit iDamineli. Conti fc Lopes! 
ll997l:lDavidsonlll99it . The 5. 52- year periodicit y was also 
found at centimeter and millimeter wavelengths lICox et all 
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[ 19951: lAbraham fe Daminelil Il999l : iDuncan fc White! 120031 : 
I Abraham et alJ l2003l) ~ VLBI observations showed a very 
compact emitting region, coincident with the star position at 
the minimum in the light curve, whi ch grew to an extended 
and e longated feature at maximum jPuncan. White fc Liml 
I 19971) . indicating the probable existence of an ionized disk 
of variable size. Infrared observations also showed a com- 
pact disk, with a power law spectrum between 2 and 8 /im, 
attributed to thermal e mission by dust at temperatures be- 
tween 1500 and 400 K jMorris et al.lll999l) . 

A draw back in the binary system hypothesis came 
from the high resolution HST spec troscopic observations 
of rj Carinae iDavidson et alJ |2000F) . which showed that, 
although the radial velocity of the emission lines formed 
close to the star varied at the different phases of the bi- 
nary orbit, they were incompatib l e with the orbital param- 
eters derived by iDamineli et alJ il997t) . Without the as- 
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sumption of a binary companion, the spectroscopic events 
could be expla ined if they were the result of periodic shell 
ejections llZanella. Wolf fc Stahllll984l:lDavidson et alll99fll 
iMartin fe Koppelmanl l2004h . Under these conditions, the 
dense gas in the shell would screen the stellar UV radia- 
tion, responsible for the formation of the ionized envelope 
surrounding r\ Carinae, favoring also the formation of dust, 
which could explain the enhanced IR emission. 

On the other hand, a strong argument in favor of the 
binary hypothesis is based in the intensity of the observed 
X-ray emission, which cannot be explained by shell ejec- 
tion but can be easily ac counted for by free- free emission 
from wind- wind collisions iPittard et al]ll998f) . In fact, this 
process was already postulated in other massive binary sys- 
tems to explain the differences between the observed X-ray 
emission (~ 10 35 erg s" 1 ) and wha t is expected from indi- 
vidual stars (~ 10 32 - 10 33 erg s" 1 ) iPrilutskii fc Usovll97ri: 
|Pollocklll987tlZhekov fc SkinneilEoOOl) . 

The 1997 and 2003 events, characterized by minima in 
the light curves at other wavelengths, were also observed 
at X-rays, preceded by a strong increase in the flux den- 
sity. The enhanced emission was attributed to an increase 
in the wind-wind collision strength as the stars approach 
each other during periastron passage; the sharp decrease to- 
wards the minima was assumed to be produced by an in- 
crease in the H column density, perhap s as the result of an 
increase in the r\ Carinae mass loss rate jPittard et alJll998t 
ICorcoran et alJl2001al : |Pittard fc Corcoranll2002l) . 

In this paper we conciliate the two hypothesis: a binary 
system and a shell ejection-like event. We use the binary 
system to explain the X-ray emission by shock formation in 
the wind-wind collision. What we call the shell ejection-like 
event is the appearence of dense material between the ob- 
server and the binary system, formed by the stellar wind 
matter accumulated behind the cool expanding shock rem- 
nant. For this process to be effective in absorbing the X- 
rays, the periastr on must be in the opo site d irection of what 
was pro posed bv lDamineli et all i200(J) and lCorcoran et alJ 
being compatible with the HST observations re- 
ported bv lDavidson et al] (120001) . We show that in this event 
dust can be formed in the shock interface, where the gas 
cools down by radiation and it is shielded from the UV stel- 
lar radiation by the large amount of accumulated gas. The 
dense gas and dust shell will eventually expand and form 
a disk in the plane of the orbit, its internal parts will be 
ionized by the UV radiation of the two stars in the binary 
system and will be seen at radio wavelengths as free-free 
emission. The disk will also slow down the wind of r\ Cari- 
nae in the plane of the orbit, which could contribute with 
other effects (e.g. non symmetrical wind from the massive 
rotating pri mary star) t o the latitude-dependent behavior 
observed bv lSmith etafl (120031) . 

In Section 2 of this paper we calculate the physical con- 
ditions in the shocked region, taking into account radiative 
absorption and emission. In Section 3 we calculate the grain 
formation rate and in Section 4 we compare our results with 
the optical and U V observations of the 2003 lo w excitation 
event reported bv lMartin fc Koppelmanl i2004h . Finnaly in 
Section 5 we present the derived X-ray light curves for dif- 
ferent values of the orbital parameters and of the r\ Carinae 
mass loss rate and in Section 6 we present our conclusions. 



2 THE SHOCK MODEL 

In massive binary stellar systems, in which both stars 
present high mass loss rates in the form of supersonic 
winds, a shocked region is formed between the two stars, 
resulting in an increase of gas density and temperature. 
One of the main observational signatures of this process is 
the strong X-ray emission, since in supersonic shocks the 
gas temperature increases from 10 4 K to ~ 10 7 — 10 8 K 
and free- free emission becomes very important; this kind 
of process is no w commonly associated to wind col lisions 
in many objects l|Corcoran et a" ] |2001atlThajier et al 1 l200ll : 
iLepine et al.ll200lh . The strength of the shock depends on 
the separation between the stars and it can increase dramat- 
ically during periastron passage in highly eccentric orbits, as 
seems to occur in the r\ Carinae system. 

For the study of the shock region we used, as in Chle- 
bowski (1989) and Usov (1992), the basic hydrodynamics 
equations for mass continuity, fluid momentum and energy 
conservation, neglecting magnetic fields. Considering an adi- 
abatic stationary shock, the jump conditions are given by: 
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where p is the gas density, v the relative wind velocity, P 
the gas pressure and 7 the adiabatic exponent; indices 1 and 
2 indicate pre and post-shock parameters. From equations 
(1) - (3) we can calculate the jump relations for density and 
temperature: 
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where M = Vi (pfkTi/ prun)" 2 ' is the Mach number, p the 
molecular weight and itlh the mass of the H atom. 

The shock front is formed at a distance ri from the 
primary star given by: 

D 



n = r (6) 

1 + n 2 

where D is the distance between the two stars and r\ = 
M a v s /MpVp. M p and M s are the mass loss rates of 77 Carinae 
and the companion star, respectively, and v p and v s their 
assimptotic wind velocities. 

For the 77 Carinae system, no direct information is avail- 
able about the individual masses, compositions and mass- 
loss rates since, even with the very high spatial resolution 
achieved by the HST, it was not possible to sepa rate the stel- 
lar spe ctra from that of the surrounding nebula. THillier et alJ 

(2001) , analyzing the P Cygni velocit y profiles, favored a 
mass- loss rate of 10 -3 Mq yr -1 , while iPittard fc Corcoranl 

(2002) found that a mass-loss rate of 2.5 x 10" 4 M yr -1 
fitted better the X-ray spectra obtained with the Chandra 
grating spectrometer. Observations indicate that r\ Carinae 
has a terminal win d velocity of v p ~ 500 — 700 km s 
dHillier et al 1 l200ll): its companion is probably an OB or 
WR type star iDamineli et alJl200ol) . with typical wind pa- 
rameters, e.g. M s ~ 10" 5 Mq yr" 1 and v s ~ 3 x 10 3 km 
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s" 1 {Abbott et. aljfl98^: iLamerslEooH: iPittard fc Corcoran! 

The orbit a l p arameters are also uncertain. While 
iDamineli et alj jl997f) obtained an eccentricity e = 0.63 
from t he radial velocity of the Pa,j and Pa<5 lines. iDavidson I 
il997f) . analyzing the same data found that an eccentricity 
of 0.8 fitted better the radial v elocity curve. An even h igher 
eccentricity was postulated bv lCorcoran et all j2001al) (e = 
0.9) to reproduce the X-ray light curve. 

Therefore, to calculate the physical parameters at the 
shock region, we will assume e = 0.9, M p = 2.5 x 
1O" 4 M yr" 1 , M s = 10" 5 M Q yr _1 , v p = 700 km s" 1 and 
v s = 3000 km s" 1 . Under these conditions, we find that at 
periastron n ~ 1.1 A.U., n\ p ~ 5 x 10 10 cm" 3 for the n 
Carinae wind and n\ a ~ 1.4 x 10 9 cm -3 for the wind of the 
secondary star, where m = pi/fimH is the gas number den- 
sity. From equations (4) and (5) and assuming Tj ~ 10 4 K 
for both winds, the post-shock density and temperature for 
r\ Carinae and the companion star winds turned out to be 



n 2p 



2 x 10 11 cm" 3 , n 2s 



5.6 x 10 a cm" 3 , T 2p 



6 x 10° 



K and T 2s ~ 2 x 10 s K. These results are compatible with 
the limits imposed to the density and tempe rature by the 
recent Chandra X-ray spec tral observations (ISeward et alJ 
l200ll:ICorcoran et alJEioOltl) . 

As a consequence of the high temperature and density 
at the shock region, the gas emits large amounts of radi- 
ation and cools down. This process is totally described by 
the first law of thermodynamics; considering also isobaric 
transformations for gas equilibrium we can write: 



3kp(t) dT(t) _ kT(t) dp(t) 
IprriH dt prriH dt 



+ V x (KVT) = T - A, 



(7) 



where K represents the thermal conduction, F and A the 
energy gain and loss functions, respectively. 

At high temperatures, the cooling functio n is dominate d 
by free-free radiation and can be expresed by (ISpitzerll97Sft : 



A ff = 2.4 x 10~ *'T*mn e 



(8) 



where rij and n e represent the number density of ions and 
electrons, respectively. 

At lower temperatures line transitions become impor- 
tant; the cooling function due to a transition j — > k of energy 
Ejk, for an element with abundance a% is: 



Aii nc „ / -, . n e 7fej 
= aiXjE jk ^ jk 1 + 



A 
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where Xj is the element degree of ionization, 7^ and •ykj are 
the collisional excitation and de-excitation parameters and 
Akj is the spontaneous emission rate. An approximate fit for 
the line cooling function was given by Matthews & Doane 
(1990): 
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where 61 = 1.53 x 10" 27 , b 2 = 1.25 x 10" 9 , p = 1.2 and 
q = 1.85. 

The photoionization h eating function (erg s" 1 cm ) is 
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Figure 1. Temperature of the shocked g function of time 

for three values of the initial gas density: 10 9 , 10 10 and 2 X 10 11 
cm -3 . In all cases the equilibrium temperature is 7500 K. 



where a v ~ 7.9 x 10" 18 Z" 2 (i/o/^) 3 <? is the photoionization 
cross section, Z is the atomic number, g is the Gaunt factor, 
vq is the Lyman limit frequency and J„ is the photon mean 
intensity. 

The time dependence of the gas temperature, deter- 
mined from Equation (7), is shown in figure Q for an initial 
temperature of 2 x 10 8 K and several values of the shock gas 
density; the equilibrium temperature of the ionized region 
turned out to be in all cases ~ 7500 K. For lower initial 
temperatures the cooling time is even shorter. 

Clearly, at this temperature dust grains cannot be 
formed or survive. However, we will show that the two stars 
cannot provide enough ultraviolet photons to ionize all the 
material that accumulates behind the shock. The neutral gas 
will then act as a shield behind which dust can be formed. 

To quantify this statement we will consider that the 
matter from the r/ Carinae wind accumulates behind the 
shock during the cooling time (about 1 day), which corre- 
sponds, at periastron, to a thickness of the accumulated gas 
of about 10 12 cm. 

On the other hand, the thickness of the ionized region 
Ar can be calculated, using typical O type stellar parame- 
ters, solving the equation: 



hv 



dv 



Awr'n^nuor' 'dr 



(12) 



where R* is the stellar radius, F v the flux density at the stel- 
lar surface and a^ 2 ' the hydrogen recombination coefficient. 
The solution of Equation (12) gives Ar ~6x 10 9 cm, much 
smaller than the thickness of the material accumulated be- 
hind the shock front. The region will have a stratified struc- 
ture, as shown in figure [5] At the borders (a), the stellar 
radiation heats and ionizes the gas. Just after this region 
(b) the high opacity maintains the gas neutral. Equation 
(7) will still give the equilibrium temperature, but with the 
photon intensity in Equation (11) given by: 



4ttJ„ = F„(R, 



Rle- T » 



(13) 



where r is the distance from the star and t„ the optical 
depth at frequency v. 
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Figure 2. Schematic view of the stratified shock region: (a) ion- 
ized region, (b) transition neutral region, (c) dust forming region 

The equilibrium temperature will go from ~ 7500 K, 
at the boundary with the ionized region, to ~ 100 K at the 
center of the shell (region "c" in fig. 2). In these conditions 
dust grains may be formed and also grow. 



3 GRAIN FORMATION 

Grain formation is not a totally understood process, it 
envolves calculations over several elements and species, 
coupled to each other. To compute the kinetic equations 
of grain composition and evolution, the nucleation the- 
ory of moments is nee ded jGail. Keller fc SedlmaviJl 19841 
iGail fc SedlmavJl987l) . However, just to estimate the mean 
particle size and the growth timescale, the classical nucle- 
ation theory is a good approach (Hoyle & Wickramasinghe 
1991). The dust growth and evaporation rates, described by 
this theory, are given by: 

(*) =*J-(M?L) (14) 
\dt) gr ^Ap s \Tipm„) ^ A h 

where a is the grain radius, p a the specific density of the 
grain material, Ak its atomic mass number and £ the sticking 
coefficient. The evaporation rate is given by: 

( — \ — ( kT dwt m <iu3t \ 5 , , 

\dt) ev ~ p s kT\ 2tt J ' 1 ' 

where P is the evaporation pressure. 

The chemical composition of the material that will form 
the grains is not known. T he nebula around r? C arinae is 
N rich and C and O poor iDavidson et al]ll986f) : infrared 
observations show the caracteristic 8-13 /im silicate feature 
iRobinson et al .119871) . However, HST observations detected 
CII and CIV lines, probably formed in the wind of the 
secondary sta r, showing a carbon-rich inner environment 
llLamers|l200ll) . For that reason, we will analyze the forma- 
tion rates of both carbonaceous and silicate grains. 

Experimental data show that carbon dust (C n ) con- 
densates at temperatures of about 2700 K at the pressure 
Psat ~ 10 4 dyn cm and that nucleation is efficient if 
P > Psat (Draine & Salpeter 1979; Gail, et al. 1984; Gail 
& Sedlmayr 1987; Andriesse, Donn & Viotti 1978; Hoyle 
& Wickramasinghe 1991). This condition is fully satisfied 
in the neutral gas region of the shell. Applying equations 
(14) and (15) for the post-shock equilibrium values we ver- 
ify that equation (14) dominates over equation (15), allow- 
ing dust growth in the internal region (part "c" of figure 
2). For large grains, evaporation becomes important and an 



equilibrium size for the dust is obtained from the condition 
yXdq/rft) = 0. For graphite, using p s ~ 2.2 g cm" 3 (D raine 
& Lee 1984), n c /n H ^ 10" 4 (ISavaee fc Sembachlll996l) and 
Ac = 12, we obtained a ~ 0.1 pm. The growth occurs in 
a time scale given by t gr ~ a/(da/dt) gr , which is of ~ 5 h. 
It means that, after the shock is formed, the time needed 
for gas cooling and further dust formation and growth to 
aproximately 10 -5 cm is of several hours. 

Considering now silicate dust formation (SiO), the satu- 
ration pressure is P s ~ 10 2 dyn/cm 2 , which gives T s ~ 1400 
K (Lefevre 1979). Considering a density p s ~ 2.5 g cm~ 3 
and using equations 14 and 15 we find that the grains form 
with a mean radius of a ~ 0.09 pm in a timescale of ~ 10 
h. These results are similar to those obtained considering 
amorphous carbon dust. 

Dust formation in the post-shock gas is also in agree- 
ment with recent observations of the binary system WR140, 
which shows periodic events at each periastron passage 
(Monnier, Tuthill & Danchi 2001). The same occurs in 
WR106, WR104 and WR137 where an increase and conse- 
quent decay in IR flux is observed llCohen & Voecl 1978; 
i Pitault et alJll983t ICohenl Il995t iKwok Volk fc Bidelman l 
1997j: IWilliams. Kideer fc van der Huchd 1200 it iKato et all 



4 SHELL-LIKE EFECTS 

The grain population formed during periastron passage must 
be the main source of visible and UV radiation absorption. 
Only recently, HST observations reached enough spatial res- 
olution to separate the stellar system from the surrounding 
nebula. During the 2003 low excitation event, the data pre- 
sented bv lMartin fc Koppelmar] i2004l) showed a wavelength 
dependent dip in the light curves for A2200, 2500, 3300 and 
5500 A, the corresponding absorption was 0.5, 0.38, 0.17 
and 0.01 mag, compatible with the interstellar absorption 
law with R = Ay /E(B - V) ~ 3 and E(B - V) = 0.08 
mag iFitzpatrickl Il999l) . Notice that the high absorption 
at 2160 A in the interstellar med ium is mainly due to 
C (iMathis. Rumpl fc Nordsiecklll977h . which is underabun- 
dant in the rj Carinae wind, b ut is probably pre sent in the 
wind of the companion star (Hil lier et al.l l200lT) . Another 
explanation for the larger decrease in the observed flux at 
the lowest wavelength, which mimics absorption by inter- 
stellar dust, is the decrease in the strength of the Fell emis- 
sion lines that also occurs during the low excitation events 
jGull. Davidson fc Ishibashl [2000 ) , in which case a lower 
amount of carbonaceous grains would be needed to explain 
the decrease in the emission. 

If the dip in the observed light curves were due to ab- 
sorption by dust with properties similar to those of the 
interstellar dust, the needed H column density would be 
( Bohli ng. Savage fc DrakeHl978T) : 

N H = 5.5xl0 21 E(B-V) em'mag" 1 = 4.6 x 10 20 cm" 2 (16) 

As we saw in Section 2, the H column density at the 
shock is about an order of magnitude higher than what is 
needed to explain the UV light curve dip, implying that ei- 
ther the fraction of dust formed at the shock is only 10% 
of the amount found in the interstellar medium, or that the 
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dust formation region covers only a fraction of the UV emit- 
ting region. 

Another consequence of the shell-like effect is the behav- 
ior of t he Ha line and its adjacent c ontinuum, as reported 
also by iMartin fc Koppelmanl <2004 . While the A6770 A 
continuum did not show any signs of absorption during the 
low excitation phase, as expected from the low optical depth 
of the grain population at this wavelength, the Ha emission 
presented a significant decrease. In the scenario presented in 
this paper, we interpret the variation in the Ha light curve 
not as a decrease in the intrinsic emission, but as an ab- 
sorption by the neutral gas in the shock front. Under these 
circumstances, the absorbed flux AFho, integrated across 
the line profile will be: 

AF Ha = !^B Ha N 2 F Ha (17) 
c 

where vhc is the frequency of the Ha line, Bho the Einstein 
coefficient for the transition and N2 the column density of 
H atoms in the energy level E2 = hv% % corresponding to the 
principal quantum number 2: 

N 2 = N H e- hU2/kT (18) 

From the observations of lMartin fc Koppelmanl |2004), 
F Ha ~ 5 x 10~ 9 erg cm -2 s" 1 and AF Ha ~ 10~ 9 erg cm -2 
s _1 , which implies that a column density N2 = 10 10 cm -2 is 
necessary to produce the observed decrease in the Ha flux. 
This is compatible with the inferred H column densities and 
temperatures of the order of 1000 K expected in the neutral 
region. 

The interpretation of the IR light cur ve is complicated 
by t he contribution of several sources dWhitelock et all 
1994). The amount of dust produced in our model is very 
small compared to the total amount of dust in the r\ Cari- 
nae nebula, formed probably after the major events of mass 
ejection. Besides, in our model the dust formation process is 
only effective during a short time (around periastron pas- 
sage) compared to the 5.5-year orbital period. However, 
since the predicted dust temperature is high, maybe it could 
be related to the em ission at short wavelen gths observed 
bv lMorris et alJ il99Sl> and lHonv et alJ (1200 ill , also explain- 
ing its anticorrelat ion with the low excita tion events seen at 
other wavelengths dWhitelock et al Jl994Tl . The presence of a 
dip in the near infrared light cu rves, coincident with the min- 
imum at other wavelengths llFeast. Whitelock fc Marand 
I 20011) . can be easily explained if part of the IR emission 
is produced by the free-free process, in the same region that 
originates the optically thick millimiter spectrum. 



5 THE 2-10 KEV LIGHT CURVE OF r, CAR 

The emitted X-rays will be absorbed by matter intercepting 
the line of sight, regardless if it is in the for m of free atoms 
and mo lecules or if it is condensed in grains. ICorcoran et alJ 
fcOOlal) were able to reproduce the X-ray light curve of the 
1997 low excitation event assuming wind-wind collisions and 
a highly eccentric orbit. In their model, periastron occurs 
when the secondary is near opposition with the Earth and 
the minimum in the light curve can be explained by the 
increase in optical depth due to the denser wind of the pri- 
mary star intercepting the line of sight. However, to fit the 



minimum duration, they had to postulate a strong increase 
in the r\ Carinae mass loss rate when the distance between 
the two stars decreases. 

In our work, we calculated the 2-10 keV X-rays emis- 
sion, integrated over all the shock surface at different orbital 
phases, using the expression given by Usov (1992): 

M 1/2 M 3/2 A/ 1 / 2 l/ 3 / 2 \ 

*** - 2.7 x 10- ^J^- + 0.6^ J , (19) 

where M p is the mass loss rate of r\ Car in 1O _4 M0 yr _1 , 
M s is the mass loss rate of the companion in 1O _6 M0 yr~\ 
V p and V a the wind speeds in units of 1000 km s _1 , and D 
is the separation between the stars in AU. Notice the D~ l 
dependence of the total flux, due to the increase in the shock 
surface size as the distance between the two stars increases. 

However, we assumed that periastron passage occurs 
when the secondary is near conjunction with the observer, 
as can be seen schematically in Figure 3. In fact, the actual 
position of the periastron relative to t he observer depends on 
how the observations are interpreted. iDamineli et al.l (Il997t) 
used the periodic radial velocity variations of the Pa7 and 
Pa<5 emission lines to determine the orbital parameters, as- 
suming that they were formed in the r\ Carinae wind. If, on 
the contrary, th e lines are formed in th e cooling, post-shock 
expanding gas llHill et alJI 2'ob"bl I 2002ft , the periastron will 
be in the opposite direction, a s assumed in our wor k and 
required by the observations of iDavidson et all |2000). 

In our scenario, as the secondary moves in its orbit, it 
leaves behind a dense and cool region, moving away from r\ 
Carinae with velocity of about 150 km s _1 . Matter from the 
primary wind will be shocked to temperatures about 3 x 10 6 
K and, after isobarically cooled, it will reach densities of 
2 x 10 16 cm -3 . Although this material will not contribute 
to the observed X-ray emission, after conjunction it will be 
positioned between the X-ray producing wind-wind shock 
and the observer, as represented by the gray tail in Figure 3. 
Notice that the absorption is due to a much larger amount of 
material than that which results from the wind-wind shock. 
The major source is the matter accumulated to the left of 
the secondary, after the star has passed, due to the primary 
wind. 

To calculate the amount of gas accumulated towards 
the line of sight at each orbital phase, we assumed that the 
observer is in the orbital plane of the binary system. Proba- 
bly, this is not true (Damineli 1997), but even so, our result 
will not be changed appreciably if the the shell is thick. We 
divided the movement of the secondary star in its orbit into 
time steps and computed the amount of gas intercepting the 
line of sight at each step, taking into account the wind contri- 
bution from 77 Carinae and the fact that the shell moves away 
from the star. The shell evolution is schematically shown in 
the lower part of figure 3, phase 1 represent the starting 
point of our calculations; before that, the shell stays always 
behind the shock, and no absorption is observed. At phase 
2, the secondary is in conjunction, and at phase 3 it reaches 
periastron; at this point part of the expanding shell already 
intercepts the line of sight and the X-rays are absorbed; at 
phases 4 and 5 the column density and X-ray absorption 
increase. Eventually, the shell expands enough to return the 
absorption to its initial value. The resulting column density, 
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Phase 
0.98 





To the Earth 



Figure 3. Evolution of the gas accumulated behind the shock, 
close to conjunction and periastron passage (left) and far from 
them (right). The upper panel represents the shocked surfaces and 
lower panel, the accumulated gas (gray tail) at different points in 
the orbit. The numbers under the phase labels represent the true 
orbital phases, calculated for an eccetricity e = 0.95 and a 5.52 
year orbital period. Orbital phase one corresponds to periastron 
passage, which we assumed to have occurred in December 16, 
1997. 
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Figure 4. Neutral gas column density of the accumulated post- 
shocked material intercepting the line of sight at different epochs 
around periastron passage for the 1997 low excitation event. The 
orbital phases shown in the upper axis were calculated with the 
orbital parameters described in figure 3. The three curves, in order 
of increasing intensity, correspond to mass loss rates for r] Carinae 
of 10" 4 , 2 X 10~ 4 and 5 X 10" 4 Mq yr" 1 . 



as a function of epoch and phase, is shown in figure 4 for 
three values of the r\ Carinar mass loss rate. 

The model X-ray light curves, taken into account ab- 
sorption, tog ether with the data for the 1997 low exci- 
tation event jlshibashi et al.l |l999T> are shown in figure 5. 
In all cases the mass loss rate of r\ Carinae was taken as 
2.5 x 10 -4 Mq yr _1 ; in the upper panel two values for the 
eccentricity were used: 0.9 (light line) and 0.95 (heavy line), 
assuming that the periastron passage coincides with con- 
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Figure 5. Model light curves for the 199 7 low excitation event , 
superimposed to the observational data Jlshibashi et al 
for a mass loss rate of 2.5 X 10 — 4 Mq yr — 1 . Upper panel: angle 
between conjunction and periastron <J> = 0°; eccentricities 0.9 
(light line) and 0.95 (heavy line). Lower panel: eccentricity e = 
0.95 and angles $ = 0° (heavy line) and <E> = 30° (light line). 



junction ($ = 0, where $ is the position angle of the sec- 
ondary in its orbit, measured from r\ Carinae); in the lower 
panel, the eccentricity was taken as 0.95 and the angle be- 
tween periastron and conjunction was varied from $ = 
(heavy line) to <E> = 30 (light line). We see that in all cases 
the model reproduces fairly well the observations. 

In the upper panel of figure 6 we compare the observa- 
tions with model light curves for three values of the mass 
loss rate: 1.0 x 10~ 4 , 2.5 x 10 -4 and 5.0 x 10~ 4 M yr _ \ 
with v p = 700 km s _1 and an orbit with e = 0.95 and 
$ = 0. In the lower panel we show the model X-ray light 
curves for two values of the primary wind velocity : 700 km 
s _1 (heavy curve) and 500 km s" 1 (light curve) for a pri- 
mary mass loss rate of 2.5 x 10 -4 Mq yr _1 and orbital 
parameters e = 0.95 and <!> — 0. As we can see, the model 
light curve depends strongly on the mass loss rate and it is 
almost independent of the chosen velocity. The best fitting 
occurs for M v = 2.5 x 10~ 4 Mq yr _1 , the value obtained by 
IPittard fc Corcoran! i2002f) based in numerical simulations. 
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Figure 6. Model light curves for the 1997 low e xcitation event, 
super imposed to the observational data data Jlshibashi et alj 
1999) for eccentricity e = 0.95, angle <J> = 0°. Upper panel: mod- 
els for v p = 700 km s — 1 and mass loss rates of 10 — 4 , 2.5 X 10 — 4 
and 5 X 10 -4 Mq yr — 1 . Lower panel: models for mass loss rate 
2.5 X 10 -4 Mq yr _1 and velocities 500 km s _1 (light line) and 
700 km s — 1 (heavy line). 



Although the calculations presented here are valid close 
to periastron passage, they should be able to reproduce the 
X-ray light curve at all phases. As the secondary moves in 
its orbit, the amount of material accumulated behind the 
shock intercepting the line of sight decreases, as shown in 
Figure 4. On the hand, as the star approaches opposition, 
the column density of material intercepting the line of sight 
due to the regular r\ Carinae wind increases. However. lUsovl 
( 1992) showed that the average hydrogen column density 
that produces the X-ray absorption decreases as the separa- 
tion between the star increases because, although absorption 
is large at the shock center, the larger emitting area more 



than compensates this effect, as can be seen in the following 
equation: 



N H ~ 1(T 



M„ 



1O- 5 M yr" 1 I \ 10 13 cm 



D 



V v 



10 3 kr: 



(20) 



In our scenario, at opposition D ~ 30 AU, resulting 
in a column density Nh ~ 10 22 cm" 2 , too low to produce 
significant absorption. 



6 CONCLUSIONS 

In this work we showed that the wind-wind collision model 
of a binary system in a highly eccentric orbit can explain 
the periodic events in the 2-10 keV X-ray light curve of r\ 
Carinae, as well as the shell-like events at optical and UV 
wavelengths, without any increase in the prim ary mass loss 
rate, as postulated bv lCorcoran et alJ (l2001al). 

As in p reviou s works Pittard et alJ Il998l : 
ICorcoran et~aTI l2001al : iPittard fc Corcoran! 12002ft . we 
calculated the physical conditions in the shock region close 
to periastron passage, for an eccentricity e = 0.9, mass loss 
rates of M p = 2.5 x 10~ 4 M yr -1 and M s = 10" 5 M c . 



km s" 1 . We obtained for the X-ray emitting region a density 
ri2s ~ 5.6 X 10 9 cm -3 and temperature T ~ 2 X 10 8 K, in 
agreement with recent Chandra observations and numerical 
wind collision models jPittard fc Corcor an 2002). 

The difference with the cited papers is that in our work 
we took into account that the wind from r\ Carinae accumu- 
lates behind the shock, forming a shell about 10 12 cm thick. 
We then showed that, as the secondary moves in its orbit, 
the shell cools due to free-free and line emission. Also, the 
number of ionizing photons from r\ Carinae and its compan- 
ion star are only enough to maintain the gas temperature at 
its equilibrium value of about 7500 K in the external parts of 
the shell (10 9 cm in depth), leaving the central part neutral 
and cold. We calculated that in these conditions dust grains 
can form and grow in timescales of a few hours. These grains 
would absorb the optical and UV radiation emitted by the 
shock material as it cools down. We used the recent HST 
observations of the 2003 low excitation event, presented by 
iMartin fc Koppelmanl (|2004) , and showed that the decrease 
in the continuum flux at wavelengths of 2200, 2500, 3300 
and 5500 A can very well be explained in this scenario. 
On the other hand, the decrease in the observed Ha emis- 
sion, without any similar effect in the adjacent continuum, 
was attributed by us to absorption by the neutral and cool 
H gas accumulated in the central region of the shell. From 
the amount of absorbed radiation we were able to calculate 
the number of hydrogen atoms in the second excited energy 
level (~ 10 10 cm" 2 ). This is compatible with the inferred H 
column densities and temperatures of the order of 1000 K 
expected in the neutral region. 

We also calculated the total amount of gas that inter- 
cepts the line of sight for each point in the orbit near pe- 
riastron passage, assuming it is close to conjunction, and 
de rived a profile similar t o that found phenomenologically 
bv lCorcoran et~ ] fcOOlal) . but without the need of an in- 
crease in mass loss rate from 77 Carinae. 

We constructed the expected 2-10 keV light curves for 
several values of the eccentricity, phase of opposition, mass 
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loss rates and wind velocities for r\ Carinae and we were able 
to fit very well the shape and duration of the emission dip. 
We found that the mass loss rate is the only critical param- 
eter that affects the shape of the X-ray light curve, since for 
high eccentricities $ angles lower than 30 degrees represent 
less than 1 day in the light curve timescale. The best fitting 
was obtained for M p = 2.5 x 1 0~ 4 yr~\ which coin- 

cides with the value obtained bv lPittard fc Corcoran! (I2002T) 
in their numerical simulations. 

An interesting result of our model is that it explains 
the difference in the duration of the dip in the light curves 
at different wavelengths. At optical and UV wavelengths the 
dip is due to absorption by dust. Since dust growth rates are 
very sensitive to gas density, its formation rate, as well as 
the absorption it causes, will decrease very fast as the dis- 
tance between the stars increases after periastron passage. 
X-rays, on the contrary, are absorbed by neutral gas, which 
will accumulate across the line of sight when the gas shell 
formed behind the shock expands away from the stars, and 
will remain there for a long time after grain formation had 
stopped. 

In the scenario presented here, the material accu- 
mulated in the shell will expand into the orbital plane 
and form a thick disk, which can be se en at millimeter 
wavelengths through i ts free - free emission dCox et alJll995t 
lAbraham fc Daminelil Il999t iDuncan fe White! I2003T) . The 
dip in the millimeter light curves can be explained by a 
decrease in the number of available ionizing photons, which 
are absorbed by the neutral shell close to the binary orbit. 

Finally we must mention that what we presented here is 
a very simplified version of a very complicated process and 
any refinement in the derived orbital and physical parame- 
ters of the binary system should be obtained from detailed 
numerical simulations. 
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